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We present a measurement of the angle fa of the CKM Unitarity Triangle using time-dependent 
Dalitz analysis of D — > K%-k + -k~ decays produced in neutral B meson decay to a neutral D meson 
and a light meson (B° —> D'*'/i°). The method allows a direct extraction of 2fa and, therefore, 
helps to resolve the ambiguity between 2 fa and n — 2 fa in the measurement of sin20i. We obtain 
sin 20i = 0.78 ± 0.44 ± 0.22 and cos 201 = l.STto.tl-oil- The si S n of cos2fa is determined to be 
positive at 98.3% C.L. 

PACS numbers: ll.30.Er, 12.15.Hh, 13.25.Hw, 14.40.Nd 



Precise determination of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix elements |l| is important to 
check the consistency of the Standard Model (SM) and 
search for new physics. The value of sin20i, where fa is 
one of the angles of the Unitarity Triangle, is now mea- 
sured with high precision: sm20i = 0.725 ± 0.037 HQ. 
This leads to four solutions in fa: 23°, 67°, (23 + 180)°, 
and (67 + 180)°. Resolution of this ambiguity has been 
attempted using time-dependent angular analysis in the 
B° — > J/ipK* a (Kgir°) decay. This technique provides a 
measurement of cos 2 fa and therefore helps to distinguish 
between the solutions at 23° and 67° 

A new technique based on the analysis of B° — > 
D[KgTr + ir~]h° has been recently suggested 0. Here we 
use h° to denote light neutral mesons, 7r°, r/ and cu. The 
neutral D meson is reconstructed in the KgTT + TT~ decay 
mode, its resonant substructure has been measured |7J,|8|. 

Consider a neutral B meson that is known to be a B Q 
at time ttag- At another time, t S i S , its state is given by 

\B°(At)) = e H At l/ 2 ^° x (1) 
(jB°) cos(AmAt/2) - i- \B°) sin(AmAt/2)) , 

where At — t s - lg — ttag, Tbo is the average lifetime of 
the B° meson, Am, p and q are parameters of B a -B° 
mixing. Here we have assumed CPT invariance and ne- 
glected terms related to the lifetime difference of neutral 
B mesons. In the SM, \q/p\ = 1 to a good approximation, 
and, in the usual phase convention, arg(p/q) = 2 fa. 

The B — > Dh° decay amplitude is dominated by the 
CKM favored b — ► cud diagram as shown in Fig. ^ with 
roughly a 2% contribution from the CKM suppressed 
b — ► ucd diagram. Ignoring the latter, a neutral D me- 
son produced in a B° decay is a D°, while that pro- 
duced in a B° decay is a D°. The D meson state pro- 
duced at time At is then given by cos(AmAi/2) - 




d < < d 

FIG. 1: Diagram for the dominant color-suppressed ampli- 
tude for B° -> Dtt°. 



«e 2i<?il ^o(-l)' |D°)sin(AmAt/2), where we use £ h o to 
denote the CP eigenvalue of h°, and I gives the orbital 
angular momentum in the Dh° system. In the case of 
B° — > D*h°, an additional factor arises due to the CP 
properties of the particle emitted in the D* decay (either 
D* -> Dtt° or D* -> D-y) Q. 



We follow Ref. |3| and describe the amplitude for a 
D° Kgir + Tr~ decay as f{rrvL,m 2 _), where m\ and 
m 2 _ are the squares of the two-body invariant masses of 
the KgTr + and Kgir~ combinations. Assuming no CP 
violation in the neutral D meson system, the amplitude 
for a D° decay is then given by f '{m 2 _ , m 2 ^) . The time- 
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dependent Dalitz plot density is defined by 



, 9 , A , e -|At|/T B o F f m 2 2 S 

P(m 2 + ,m 2 _,At,q B ) = — +' "M l + gg x 

oTgo ZiV V 



2 2 



) cos( AmAt) + S(m 2 _ , to^ )s'm(AmAt)} 



A= (|/(m 2 _,m 2 )| 2 - \f(ml,m 2 _)f)/F(ml,m 2 _), 
-2£ h0 (-l) l Im{f(m 2 _,m\)f* (ml , m 2 _ ) e 2 ^ 1 } 
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F(ml, m?_) 
F = \f(m 2 _,m 2 + )\ 2 + \f(m 2 + ,m 2 _)\ 2 , 

N= I \f(m 2 _,m 2 + )\ 2 dm 2 + dm 2 _, 



(2) 



where the 6- flavor charge is qs = +1 (— 1) when the 
tagging B meson is a B a (B°). Thus the phase 2<f)i can 
be extracted from a time-dependent Dalitz plot fit to 
B° and B° data if f(rn 2 +1 m 2 _) is known. Note that this 
formulation assumes that there is no direct CP violation 
in the B decay amplitudes. 

This analysis is based on 386 x 10 6 BB events collected 
with the Belle detector at the asymmetric energy e + e~ 
collide r lid . The Belle detector has been described else- 
where |llj. We reconstruct the decays B° — > Dh° for 
h° = 7T°, rj and u and B° -> D*h° for h° = tt° and 77. 

Charged tracks are selected based on the number of 
hits and impact parameter relative to the interaction 
point (IP). To reduce combinatorial background, a trans- 
verse momentum of at least 0.1 GeV/c is required of each 
track. All charged tracks that are not positively identi- 
fied as electrons are treated as pions. 

Neutral kaons are reconstructed via the decay K s — > 
7r + 7r~. The 7T7T invariant mass is required to be within 
9 MeV/c 2 (~ 3<t) of the K° mass, and the displacement 
of the 7T + 7r - vertex from the IP in the transverse (r-tp) 
plane is required to have a magnitude between 0.2 cm 
and 20 cm and a direction that agrees within 0.2 radians 
with the combined momentum of the two pions. 

Photon candidates are selected from calorimeter show- 
ers not associated with charged tracks. An energy de- 
position of at least 50 MeV and a photon-like shape are 
required for each candidate. A pair of photons with an 
invariant mass within 12 MeV/c 2 (2.5<r) of the 7r° mass 
is considered as a 7r° candidate. 

We reconstruct neutral D mesons in the K s tt + tt~ de- 
cay channel and require the invariant mass to be within 
15 MeV/c 2 (2.5a) of the nominal D° mass. D*° can- 
didates are reconstructed in the D°7r° decay channel. 
The mass difference between D*° and D° candidates is 
required to be within 3 MeV/c 2 of the expected value 
(3a). u> candidates are reconstructed in the 7r + 7r~7r° 
decay channel. Their invariant mass is required to be 
within 20 MeV/c 2 (2.5 T) of the uj mass. We define the 
angle 9 U between the normal to the uj decay plane and 
opposite of the B direction in the rest frame of lu and 
require |cos# w | > 0.3. We reconstruct 77 candidates in 



the 77 and 7r + 7r~7r° final states and require the invari- 
ant mass to be within 10 and 30 MeV/c 2 (2.5a) of the 
77 mass, respectively. The photon energy threshold for 
the prompt ir° and 77 candidates coming from B decays 
is increased to 200 MeV in order to reduce combinato- 
rial background. We remove 77 candidates if either of the 
daughter photons can be combined with any other pho- 
ton with E 1 > 100 MeV to form a 7r° candidate. 

We combine either D and h a = {ir°, 07,77} or 
D* and h° = {7r°,?7} to form B mesons. Sig- 
nal candidates are identified by their energy dif- 
ference in the center-of-mass system of the T(4S) 
(CM), AE = (J2, Ei) - E heam and the beam -energy 
constrained mass, M bc = v^bcam ~ (X^) 2 > where 
-E/beam is the beam energy and pi and Ei are the mo- 
menta and energies of the decay products of the B meson 
in the CM frame. The masses of n°, 77 and candi- 
dates are constrained to their nominal values to improve 
AE resolution. We select events with Mb c > 5.2 GeV/c 2 
and lAi^l < 0.3 GeV, and define the signal region to 
be 5.272 GeV/c 2 < M bc < 5.287 GeV/c 2 , -0.1 GeV < 
AE < 0.06 GeV (tt°, 77 -> 77) or \AE\ < 0.03 GeV (uj, 
77 — > 7r + 7T~7r°). In cases with more than one candidate 
in an event, the one with D and h° masses closest to the 
nominal values is chosen. 

To suppress the large combinatorial background domi- 
nated by the two-jet-like e + e~ — ► qq continuum process, 
variables that characterize the event topology are used. 
We require |cos$thr| < 0.80, where #thr is the angle be- 
tween the thrust axis of the B candidate and that of the 
rest of the event. This requirement eliminates 77% of the 
continuum background and retains 78% of the signal. We 
also construct a Fisher discriminant, T, which is based 
on the production angle of the B candidate, the angle 
of the B candidate thrust axis with respect to the beam 
axis, and nine parameters that characterize the momen- 
tum flow in the event relative to the B candidate thrust 
axis in the CM frame |l2| . We impose a requirement on T 
that rejects 67% of the remaining continuum background 
and retains 83% of the signal. 

Signal yields and background levels are determined by 
fitting distributions in AE for candidates in the Mb c sig- 
nal region. For each mode, the AE distribution is fitted 
with an asymmetric Gaussian for signal and a linear func- 
tion for background. The signal shape is fixed, based on 
MC simulation. The region AE < —0.1 GeV is excluded 
from the fit to avoid contributions from other B decays. 
The results from our fits to the data are shown in Figure|5] 
and Table UJ We study the systematic error of the fit by 
varying the shapes for signal and background and chang- 
ing the fit range. The difference in the signal yields does 
not exceed 5%. We also confirm that there are no feed 
across between channels and other peaking background 
by using generic BB MC. 

The signal B decay vertex is reconstructed using the D 
trajectory and the IP constraint. The tagging B vertex is 
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TABLE I: Number of events in the signal region (N to t), de- 
tection efficiency, number of signal events from the AE fit 
(iVsig) and signal purity for the B — > D^*'h° final states. 



Process 


iVtot 


Efficiency (%) 




Purity 


Dtv° 


265 


8.7 


157 ± 24 


59% 


Duj 


88 


4.1 


67 ± 10 


76% 


Dr) 


fOf 


3.9 


58 ±13 


57% 


D*n°,D*ri 


67 




43 ± 12 


64% 


Sum 


521 




325 ± 31 


62% 



obtained with well-reconstructed tracks not assigned to 
the signal B candidate and the IP constraint 13]. The 
time difference between signal and tagging B candidates 
is calculated using At = Az/^[3c and Az = zqp — z tag . 
The proper-time interval resolution function R s i g (At) is 
formed by convolving four components: the detector res- 
olutions for zcp and Ztag, the shift in the z tag vertex po- 
sition due to secondary tracks originating from charmed 
particle decays, and the kinematic approximation that B 
mesons are at rest in the CM frame |l3j . A small compo- 
nent of broad outliers in the Az distribution, caused by 
misreconstruction, is represented by a Gaussian function. 
Charged leptons, pions, kaons, and A baryons that are 
not associated with a reconstructed B° — ► D[Kgir + ir~]h° 
decay are used to identify the 6-flavor of the accompany- 
ing B meson. The tagging algorithm is described in detail 
elsewhere ITil . 

We perform an unbinned time-dependent Dalitz plot 
fit. The negative logarithm of the unbinned likelihood 
function is minimized: 

n 

-2 \ogL = -2^2 log{(l - /bg) ^ig + /bg Pbg}, (3) 

i=l 

where n is the number of events. The function 
P s i g (m+, mi, At) is the time-dependent Dalitz plot den- 
sity for the signal events, which is calculated according to 
Eq. J5J and incorporates reconstruction efficiency, flavor- 
tagging efficiency, wrong tagging probability and At res- 
olution. The function Pb g is the probability density func- 



tion (PDF) for the background. Both P s j g and Pb g are 
normalized by / P s i g ,bg(w+, mL, At)dm 2 i _dm, 2 _dAt = 1. 
The event-by-event background fraction /b g (A£, M bc ) 
is based on signal and background levels found by fit- 
ting AE as described above, the AE shape used in the 
fit, and an M^ c shape that is the sum of a Gaussian 
signal and an empirical background function with kine- 
matic threshold and shape parameters determined from 
off-resonance data. 

We describe the background by the sum of four com- 
ponents: B decays containing a) real D mesons and b) 
combinatorial D mesons, and qq events containing c) real 
D mesons and d) combinatorial D mesons. The Dalitz 
plot is described by the function f{rn/L,m 2 _) for a) and 
c). For b) and d) we use an empirical background func- 
tion which includes enhancements near the edges of the 
Dalitz plot as well as an incoherently added K*(&%2) con- 
tribution Q . The shape of this function is obtained from 
an analysis of events in the D mass sideband. The At 
distribution for the B decay backgrounds is described by 
an exponential convolved with the detector resolution. 
For the qq background, a triple Gaussian form is used, 
which is obtained from events with | cos#thr| > 0.8. The 
use of this sideband region has been validated using MC. 
We use the experimental data and generic MC to fix the 
fractions of background components. Figure shows the 
Dalitz plot distributions for candidates in signal and M bc 
sideband region, integrated over the entire At range and 
B° and B° combined. We can see clear differences in 
these distributions. 

The procedure for the At fit is tested by extracting t b + 
using B + — > D°[KgTT + 7r^]'K + decay. We obtain t b + = 
1.678 ± 0.043 ps (statistical error only), consistent with 
the PDG 2] value 1.638 ± 0.011 ps. 

We perform a fit by fixing t b o and Am at the PDG 
values with a fixed background shape and using sin2^i, 
cos 2(f>i as fitting parameters. The results are given in Ta- 
ble [n] for each of the three final states separately and for 
the simultaneous fit over all modes. We check goodness- 
of-fit using one dimensional projections to K^ir^ and 
7r + 7r~ invariant masses and St and find no pathological 
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FIG. 3: Dalitz plot distribution for the Dh° candidates from 
B signal region (left) and Mb c sideband. 



TABLE II: 
only. 



Fit results for the data. Errors are statistical 



Final state 


sin 20i 


cos 2^1 


Dtt°, 077(77] 
Dlj, Drj[3n] 
D*Ti ,D*r 1 
Simultaneous fit 


0.43 ± 0.90 
1.07 ±1.14 
0.78 ± 0.44 


2 D7 +a78 

z - u '-0.91 

-, cq+0.67 
J--OO_0.93 



behaviour. To illustrate, the raw CP asymmetry dis- 
tribution for D^*'h° candidates with an additional con- 
straint |M-+-- — 0.77| < 0.15 GeV/c 2 , to select events 
consistent with D — > Kgp, is displayed in Fig. 0] For 
D*h° candidates we take into account the opposite CP 
asymmetry. In this case the system behaves approxi- 
mately as a CP eigenstate, with an asymmetry propor- 
tional to — sin 20i. 

Uncertainty of the D — » -ftTg7r + 7r~ decay model is one 
of the main sources of systematic error for our analy- 
sis. We repeat the fit using two additional decay models 




FIG. 4: Raw asymmetry distribution for the D w [Kgp°]h° 
candidates. The smooth curve is the result of the fit to the 
full Dalitz plot. 



from CLEO Q and similar Belle analysis [l5|. The dif- 
ference between these models and our primary model Q 
is in describing of wide resonances in 7r + 7r~ and K^n, 
non-resonant part and doubly Cabibbo suppressed chan- 
nels. The CLEO does not include wide resonances 
ct(600) and /o(1370), and the doubly Cabibbo suppressed 
channel D° -> K*+(U30)ir-. Another Belle model [3 
has an additional contributions from i4T*(1410) ± pi =F and 
450). The difference in fitted values for sin20i 
and cos 2<j>i between the nominal model || and others is 
found not to exceed 0.1, and we assign this value as a 
model uncertainty. 

We vary the background descriptions to estimate the 
systematic uncertainty due to the background parameter- 
ization. We use only a combinatorial and only a signal 
D PDF for the Dalitz plot distribution. For the time 
dependence, we consider cases with only a qq component 
or only a BB component. The differences do not exceed 
0.2 and we take this value as a systematic error. 

Other contributions to the systematic error are found 
to be small: vertexing and flavor tagging (0.02), neglect- 
ing suppressed amplitudes (0.01), signal yield determina- 
tion (0.02). 

The measurement of sin20i = 0.78 ± 0.44 ± 0.22 
is consistent with the high statistics measurement in 
the J/ipK° channel 0. The result of cos20i = 
1.87to'53_o'32 allows one to distinguish between two so- 
lutions in 0i : 23° and 67°. We define the confidence level 
at which the 67° solution (negative value of cos 20i ) can 
be excluded as CL(x) = /+(_■)/(/+ (x) + /_(_•)), where 
f+(x) (/_(«)) is the likelihood to obtain the fit result 
cos 20i = x when true cos20i value of 0.689 (—0.689). 
To evaluate /+ and /_ we use sample of 2500 pseudo- 
experiments with the same size as data for both hypothe- 
ses. We fit these distributions with a sum of two Gaus- 
sians. We calculate CL for x — 1.87, 1.55 and 2.09 to 
take into account systematic uncertainties of _n' 32 in our 
cos 20i measurement. As a final result we use the small- 
est value CZ(1.55) = 98.5 ± 0.2%, excluding the 67° so- 
lution at 98.3% C.L. 

In summary, we have presented a new method to 
measure the Unitarity Triangle angle 0i using a time- 
dependent amplitude analysis of the D — > KgTr + ir~ de- 
cay produced in the processes B° — > D^*'h°. We find 
sin 20i = 0.78 ± 0.44 ± 0.22 and cos 20i = 1.87±g;fg±g;|l. 
The sign of cos20i is determined to be positive at 98.3% 
C.L., strongly favoring the 0i = 23° solution. 
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